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Abstract 

We model the behavior of a big (Brazil) nut in a medium of smaller 
nuts with a stochastic asymmetric simple exclusion dynamics of a polymer- 
monomer lattice system. The polymer or 'rod' can move up or down in 
an external negative field, occupying N horizontal lattice sites where the 
monomers cannot enter. The monomers (at most one per site) or 'fluid 
particles' are moving symmetrically in the horizontal plane and asym- 
metrically in the vertical direction, also with a negative field. For a fixed 
position of the rod, this lattice fluid is in equilibrium with a vertical height 
profile reversible for the monomers' motion. Upon 'shaking' (speeding up 
the monomers) the motion of the 'rod' dynamically decouples from that 
of the monomers resulting in a reversible random walk for the rod around 
an average height proportional to log N. 
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1 Introduction. 



Studying the coupled dynamics of granular matter of different shapes and sizes is 
of great interest for a range of phenomena. One example is the size segregation of 
particles as a result of vibrations. A typical realization is a can with nuts; upon 
shaking the larger (often taken to be Brazil) nuts rise to the top. Because of its 
wide interest the phenomenon has been considered and reconsidered and while 
some of the aspects are well-understood not everything has stopped surprising. 
If one asks for an analysis starting at the microscopic level the situation is 
not so satisfactory and even simple models have escaped serious mathematical 
handling (cf. Q, Q for further references). 

In this paper we consider such a microscopic — albeit stochastic — dynamics 
for the motion of a large particle or rod in a lattice fluid composed of monomers. 
The problem of the present paper is however not quite similar to the canonical 
Brazil nuts scenario as we are interested in the equilibrium dynamics. In fact, 
as we will see, on the time scale of the motion of the rod, the monomers are in 
equilibrium for a reversible density profile. The rod then finds its hydrodynamic 
equilibrium at a vertical height where the density of the fluid is about equal to 
its own density. Going beyond equilibrium conditions, e.g. starting from a 
homogeneous density for the lattice fluid, gives rise to additional mathematical 
problems that we will only be touching at the end of the paper (see Section [|, 
Remark 2), and which will be the subject of future work. 

The result of this paper can be classified under the heading: how to obtain a 
Markovian reduced dynamics? This problem is of course a very common one in 
noncquilibrium statistical mechanics where one considers the system composed 
of various types of degrees of freedom. The dynamics is globally defined in 
which the various degrees of freedom are coupled. In some circumstances and 
under some limit procedures one then expects that some degrees of freedom of 
the system effectively decouple giving rise to an autonomous (in many cases, 
Markovian) dynamics for a subset of degrees of freedom. In our case, it is the 
shaking, the speeding up of the monomer dynamics in the horizontal direction, 
that does the job. In this way, between any two moves of the polymer, the 
monomer configuration has the time to relax to its reversible measure and the 
polymer always sees the fluid in equilibrium. 

Our main result is a mathematically rigorous proof of this dynamical decou- 
pling between the motion of the rod and the monomer fluid when the monomer 
dynamics is (infinitely) speeded up (at least) in the horizontal direction (or- 
thogonal to the motion of the polymer). In that limit of excessive horizontal 
shaking the reduced dynamics of the polymer becomes that of a random walker 
with rates directly given in terms of the equilibrium fluid density. When N (the 
length of the polymer) is sufficiently big (depending on the rates for jumping up 
or down) the polymer finds its most probable height around its mean position 
of order log TV with a variance of order 1. In the next section we describe the 
model and the result. The third section is devoted to the proofs. The final 
section contains an open problem and some additional remarks. 
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2 Model and Results. 



2.1 Model. 

2.1.1 Configuration. 

For convenience we put the system on the square lattice Z 2 . A point i — (x,y) 
of the lattice has a 'vertical' coordinate y and a 'horizontal' coordinate x. We 
also write i — (11,12) if, in the notation of the coordinates, we want to remember 
the site i. 

The system contains a rigid polymer (large particle, rod) whose position at 
time t is denoted by Y t . For simplicity we allow the rod to move only vertically. 
The horizontal coordinate is fixed (at 0) and Y t takes values in 7L (thought of 
as the 'vertical' axis). The same results would hold if the polymer also jumps 
horizontally at rate 1. The polymer occupies N G {2, 3, . . .} lattice sites. If the 
polymer has position Y t = y, then it occupies the region 

A N (y) = {(0,y),(l,y),...,(N-l,y)}. 

This region is forbidden for the monomers (fluid particles). The monomer con- 
figuration is denoted by 77 £ {0, 1}^ and we use r\t to denote the random field 
of monomers at time t. We have that rj t (i) = if there is no monomer at site i 
at time t; i] t (i) = 1 if there is a monomer at site i at time t. The dynamics will 
always be subject to the restriction (exclusion) that rj t (i) — for i 6 AN(Y t ) 
(the rod acts as an obstacle for the fluid motion). The full configuration space 
is denoted by fl = {0, l} 2 ' x Z. 



2.1.2 Dynamics. 

We now define the coupled dynamics for the polymer-monomers system. All 
motion is via jumping to vacant sites. There are the horizontal jumps of the 
monomers (which we take symmetric and at rate 71), the vertical jumps of the 
monomers (asymmetric at rate 72) and the vertical jumps of the rod (asymmetric 
at rate 1). The asymmetry in the vertical direction models the presence of an 
external (e.g. gravitational) field acting on fluid matter and polymer but can 
in general be different for monomers and polymer. Increasing the rates 71 and 
72 can be used to simulate the greater mobility of the smaller particles upon 
shaking. We are most interested in the case where 72 — 1 and 71 ^> 1 (horizontal 
shaking) . 
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Fig.l: A polymer between monomers 
7 



Fig. 2: Rates of jumping for monomer and polymer where 72 = 7 and 71 = 1. 

Here comes the formal definition of the generators of these motions. Let / 
be a local function on £1 (i.e., a function that depends on the configuration 
in a finite region of Z 2 ). The first part of the generator represents horizontal 
monomcr-j umping : 

^(v,y) = l E I [(ij)^A N ( y ) = (D][f(v i ' j ,y)-f(v,y)] (2.1) 

(»j):*2=j2 

where rf ,3 (k) = rj(i) if fc = j, = if = i and = 7j(fc) otherwise; the sum- 
mation is over nearest neighbor pairs (ij) with the same vertical coordinate 
(*2 = 32) ■ The notation /[•] will always stand for the indicator function. 
Second comes the vertical monomer-jumping: with p < q, 

£vf(v,y) = *2 + l))/[(ii,*2 + l) i A N (y)] 

i 

x [/(»7 i>(illia+1) , »)-/(»?, y)] 

+ qrj(i) (1 - *7(*i,i 2 - 1)) J[(ii,»2 - 1) ^ ^iv(y)] 

xtffa*^- 1 )^)-/^)]}. (2.2) 

Finally, there is the polymer-jumping: with a,b £ 1R + , 

^poiy/(r?,y) = o /[»?(*) = 0,VieA JV (i, + l)][/( )7 ,y+l)-/(»7,y)] 

+ &Jfo(i) = 0,Vz e ^(y - 1)] [/(rj,y - 1) - /(»?, y)](2.3) 
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We will then choose p/q,a/b < 1 to represent an external field in the vertical 
direction driving all particles, big and small, downward. E.g. in the case of a 
gravitational field, we could have p/q = exp(— mg/kT), a/b = exp(— M g / kT) , 
where m, M denotes the mass of a monomer, resp. polymer. 
The formal generator £ of the full dynamics consists of three pieces: 

C = 7i£ h + 72>C V + £ P oiy, (2.4) 

where 71 , 72 > are additional parameters governing the rates of the monomer- 
jumping. Notice that £ mono = 7i£h + 72<Cv works on the configuration of 
monomers only (for fixed rod position), while £ po iy works on the polymer con- 
figuration (for fixed monomers). The only interaction is by excluded volume. 



The generator (2.4) can be rewritten in the form: 

Cf(r i ,y)=£ r ,f r ,(y)+£ y f y (r,) (2.5) 
where /„(•) = /(?/, ■) and /„(•) = /(•, y). 



2.1.3 Initial distribution and extra remarks. 

At time (starting time) we put the polymer at the origin: Y t= o — 0. Then fix 
a real parameter k and distribute the monomers independently with density 

K(p/qY 2 

Pil) = l + «(p/ g )* (2 ' 6) 

varying in the vertical direction (constant in the horizontal direction), condi- 
tioned on 17 (i) = 0,Vi € Ajv(0). 

More precisely, we let v p denote the product measure on {0, 1} with density 

v p (v(i)) = P(i), (2-7) 
defined by (|2.6|) . This measure is reversible for each of the mono mer ge ner ators 



process without polymer — i.e. the generators defined by (2.1) and (2.2) but 
without the indicator functions prohibiting jumps. The proof of this fact is a 
simple computation. The one-dimensional analogue is well known, see 

For any given y £ Z, we write 

v y p =MMi) = 0,Vi&A N (y)), (2.8) 
At time 0, we put the distribution (Iq on defined by 

Mo(*7> y) = $y,o v v p(dri), ( 2 -9) 
where 5 y fi stands for the Kronecker-delta. 

From the initial condition described above and the dynamics defined via 



(2.4) the process {rj tl Y t ) is generated. The measure at time t > is denoted by 
fit- Of course this depends on the choice of parameters p, q, a, 6,71 and 72 and 
we will sometimes make this explicit in the notation. 

A useful way to imagine the process is by associating two exponential clocks 
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(at rate a respectively b) to the polymer: one clock gives rise to the trial times 
for the polymer to jump up, the other indicates the trials for the polymer to 
jump down. If, just before the trial time r, say for jumping up, there are 
no monomers right above the polymer (rj T -(i) — 0,Vi G An(Y t - + 1)), then 
the jump is performed and at time r the polymer is at height Y T = Y T - + 1, 
otherwise it stays where it was. Between the trial times of the polymer, only 
the monomers move. The dynamics for the monomers for a fixed position of the 
polymer (say at y) is generated by 

^ono/W=7l^ + 72^ (2.10) 

which can be read off from ( |2.l| ) and ( |2.2| ) . The associated semigroup is denoted 
by S v (t). Now, the important thing where the 'equilibrium' in the title of this 
paper refers to, is that v v p is a reversible measure for S y (t) . This will be proven 
as Lemma [O] in Section ||. 



2.2 Results. 

2.2.1 Limiting random walk. 

In the limit 71 f +00 the motion of the rod will decouple from the monomer 
dynamics. It will be a random walk. We first introduce this limiting rod motion. 
For a, b 6 H + consider the continuous time random walk on TL with generator 

£ RW f(y) = a[l - p(y + l)] N [f(y + 1) - f(y)] + b[l - p(y l)] N [f(y - 1) - f(y)], 

(2.11) 



where the density profile p is obtained from (|2.6|). Remark that, in the notation 



of (2.5) 



£ RW f(y) = j v y P (dvKf(y) (2.12) 

£RW g enera tes a continuous time random walk Y" 4 RW which we start at = 
and with rate for moving one step upward a[l — p(y+ 1)) N and rate moving one 
step downward equal to 6[1 — p(y — 1)) N . We fix the initial state to be for 
the sake of definiteness. Our results hold for any other initial (deterministic or 
random) state. 



Proposition 2.1 Ifa/b > (p/q) N , then the random walk with generator (2.11) 
defined above has a unique reversible probability measure m on TL, which is given 
by 

m{y)= ' w (2 . 13) 

Z (l + (p/q)v) N ' 

where Z is a normalizing constant. In particular, the random walk is positive 
recurrent. 

Proof: Reversibility of m(y) is immediate, and the condition a/b > (p/q) N 
guarantees that m(y) can be normalized (i.e. Z < 00). Positive recurrence 
follows immediately from the existence of a reversible probability measure. ■ 



G 



Remark that the condition a/b > (p/q) N in the case of a gravitational field 
just means M/N < m, i.e. the density of the polymer is smaller than the density 
of the monomer-fluid. It is thus very natural that in this case the polymer will 
drift up and will float at a height where the fluid density is proportional to 1/N, 
see ( p0| ) and §. 



In order to study some global properties of the limiting random walk, in 
particular its behavior for large N, we replace the discrete distribution m(y) on 
Z by a continuous distribution: 



exp(— ax) 1 
(1 + cxp(-0x)) N Z(a,0, N) ' 



m(dx) := „ ,„^,'„,_ nT ^ x. (2.14) 



Here 



7 , , v , , exp(-ax) 1 T(a/f3)T(N - a/0) 
Z(a, 0,N) = I dx = - p m , (2.15) 



— oc 



and e a = a/b, e 13 = p/q. From (2.15) we can calculate the cumulants of the 



continuous distribution m(dx): in particular 



xm(dx) = (2.1<>) 



where ip(x) — V (x) /T(x). Using the asymtotic expansion 



m = lo gZ -j- z -^- 2+ ... (2.17) 



we obtain 

1 



xm(dx) = - logiV+ 0(1), as N -> oo (2.18) 

P 

and all higher order cumulants are of order 1 as N tends to infinity. The modus 
of m (the position where m{x) reaches its maximum) is 



Mo(m) = {-logip/q))- 1 \og{^^ 



1 



4 log TV, as N -> oo. (2.19) 



2.2.2 Main result. 

Our main result states 

Theorem 2.1 Let < p < q < oo and a, b G M + and consider the joint 



monomer-polymer process with generator (2.4)- For any finite time-interval K , 
the marginal law of the polymer motion (Y/ 71 : t £ K) converges, as 71 — » 00, 
to the law of the random walk (i^ t RW : t G K) defined by ( 2.l\ ). 
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2.2.3 Discussion. 



Since for a/b > (p/q) N the limiting motion is an ergodic random walk in a 
countable state space, the process starting from any initial distribution will 
converge to the (unique) invariant measure. Hence, by ( 23a ), the polymer will 



rise from the zero level to a level at height proportional to 4 log AT. If it starts 
in equilibrium, then it will perform a random walk around this position. This 
is exactly what we would expect from general hydr ody namics, see [||. After all, 
the fluid density at height ^ log N is precisely, cf : 



p(-\ogN) = K,N , ~- (2.20) 
confirming Archimedes' law in this model of granular matter. 



3 Proofs. 



3.1 Outline of proof 

In this section we state the main steps of the proof of Theorem 2.1. The reader 
may use this section as a guideline to the next section. The main idea of the 
proof is that in the limit 71 t 00 the monomers are moving very fast in the 
horizontal directions and thus can reach equilibrium in the time between two 
successive jumps of the polymer. Therefore the rate at which the polymer 
jumps, which is a function of the whole monomer configuration, can be replaced 
by the expectation of that rate in the equilibrium distribution of the monomer 
configuration. 

As a first step (Lemma 3.1) we identify the reversible equilibrium measure 
for the monomers for fixed position of the polymer. This is (by reversibility) the 
original reversible measure of the monomer gas without polymer, conditioned 
on having no monomers on the lattice sites occupied by the fixed polymer. 

In a second step (Lemma 3.2-Proposition 3.1) we prove that in the limit 
71 t 00 any time dependent expectation of a function f(Y t , r]t) of both polymer 
position Y t and monomer gas configuration rj t can be replaced by the expectation 
of a new function depending only on the polymer position, and obtained from 
/ by integrating out the 77 variables over the equilibrium measure. The main 
ingredients in the proof of that statement are 

1. Discrepancies in the asymmetric exclusion process move as "second class 
particles" which are a kind of random walkers. When 71 f 00, this "ran- 
dom walker" diffuses away very quickly. 

2. The distribution of the monomers at any jumping time of the polymer 
is absolutely continuous w.r.t. the monomer-fixed polymer equilibrium 
measure. 
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In the first two steps we obtain convergence of the distribution of the polymer 
position V t 71 to the distribution of the random walk Y t . To finish our proof, we 
still have to prove that the whole process {Y^ 1 : t > 0} converges to the whole 
process {Y t : t > 0} (i.e. the distributions on trajectories converge). This final 
step is made by first proving that any limiting process is Markovian and next 
that there exists a limiting process (tightness). 



3.2 Proof of Theorem 2.1 

We start this section with an easy lemma on reversible Markov processes. 

Lemma 3.1 Let {rj t : t > 0} be a Markov process on £1 with generator C and 
let \i be a reversible measure for C. Suppose A C such that fi(A) > and such 
that 1a is in the domain of the generator. Consider the process with generator 

C A f = U£(W) - (UCl A )f (3.21) 

That is, Ca corresponds to a process with "forbidden region" A c (i.e., jumps 
from A to A c are suppressed). Then the measure \iA '■= m('|-^-) * s reversible for 
Ca- 



Proof: Because the second term in the right hand side of (3.21) is just mul- 
tiplication with the function \aC(1a), it suffices to show that CAf '■= 1a£(1a/) 
defines a symmetric operator on L 2 (A, [ia)- Let /, g be in the domain of Ca- 
Since d^A — (1/ 'h(A))1a d/x, we get, using the symmetry of C in L 2 (/i): 

J g(c~Af)dnA = -j^j J iAgC A fdn 



l 



C{gl A )lAf dji 



C A gfd^ A - (3.22) 



Note that reversibility is crucial in the proof of this lemma. Indeed if fi is only 
stationary, then we cannot conclude in general that fiA will be stationary for 
the process with generator Ca- Indeed, one easily computes 

J C A fd^ A = J (UC*1 A - l A Cl A )fdfi, (3.23) 

i.e., ha will be stationary iff 1 a C*1a — IaCIa = /i-a.s. Since the profile mea- 
sures are reversible for the exclusion process of the monomers without polymer, 



we can apply lemma |j| for /i = v p , A = {-q e {0, 1}^ : J2 z eA( v ,N) v( z ) 0}, 
i.e., those monomer configurations which are excluded when the polymer is at 
vertical position y. This yields: 

Corollary 3.1 For fixed polymer position at y £ TL the measure v v p is reversible 
for the monomer dynamics with semigroup S v (t). 
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n?2 

Lemma 3.2 Fix y e Z. Let f be a local function on {0,1}^ which only 
depends on the monomer configuration in the layers at height y + 1 and y — 1 . 
Suppose that 

«£(/) = 0. 

TTien, /or any t > 0, 

lim = 0. 

7it+oo 



Proof: Abbreviate (i := and consider the case f y (rf) = 1a(t])~ MA) for a set 
A in the space of configurations depending only on a finite number of coordinates 
in labels y — 1 and y + 1. The extension to general local / is straightforward. 
Denote 

D A := {x 6 I 2 : l A (r]) + \ A (rf) for some r,}, (3.24) 
the dependence set of A. By reversibility: 

j(Sy(t)fy) 2 dH = j(SV(2t)fy)fydli. 

= M^)(^ ( . |A) (U(^))-]E^(U(^))). (3.25) 

To compute the difference of the expectations in the above expression we realize 
the processes with initial configurations r\ and ( in the same probability space 
(coupling) . 

To construct this coupling we first associate two Poisson clocks to each site 
of % with parameters 72P and 729 respectively and use them to decide the times 
of the vertical attempted jumps. A jump from (ji, 12) to (ii, 12 + 1) is performed 
at time t if an event of the Poisson process of rate p occurs at that time, a 
particle is present at (21,^2) and no particle is present at (ii,Z2 + 1) at time t— . 
Similarly, a jump from (ii , 12) to (ii, 12 — 1) is performed at time t if an event of 
the Poisson process of rate q occurs at that time, a particle is present at (11,12) 
and no particle is present at (i\,i2 — 1) at time t— . Jumps cither to or from 
sites occupied by the rod are suppressed. This takes care of the vertical jumps. 
See Ferrari (1992) for details of this construction. For the horizontal jumps 
we associate Poisson clocks with rate 71 to pairs of horizontal nearest-neighbor 
sites. When the clock associated with sites (11,12) and (i\ + 1,12) rings, the 
contents of those sites are interchanged. Also here, if at least one of the sites 
is occupied by the rod, the jump is suppressed. The horizontal motion is also 
called stirring process. See Arratia (1986) for details of this construction. More 
rigourosly, let (N t (i,j) : i = (11,12) € 1,j = (h + M2)), (N t + (i) : i £ Z) and 
(N^(i) : i <E Z) three independent families of independent Poisson processes 
of rates 71/2, pj 2 and qj 2 , respectively — the Poisson clocks. Use the notation 
dA^(-) = 1 if there is an event of the Poisson process (N t (-)) at time t, otherwise 
it is zero. The motion is defined by 

df(vt) = £ dN t (i,j)i[{ij)nA N (y) = 9}[f(v id ,y)-f(v,y)} 

(ij) -12=32 

+ £{dJV+(i) rt(i) (1 - V (k,i2 + 1)) I[(h,h + A N (y)\ 

i 

x[fW^ +1 \y)-f(n,y)] 
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+ dN t -(i) V (i) (l-7](i u i 2 - 1)) I[{k,i 2 - 1) i A N {y)} 

X [f(r 1 i ^- 1 \y)-f(r 1> y)}}. (3.26) 



Standard arguments, see for instance Durrett (1993) show that (3.26) defines a 
process t]t — ^(r/o; N[0,t]), with initial configuration 770, where $ is the func- 
tion induced by fT26| ) and N[0,t] := (N a (-,-), N+(-), N~(-) : < s < i); 



furthermore it is immediate to see that rj t has generator C y . Given two initial 
configurations rj and £, the coupling of their evolutions is constructed using the 
same Poisson processes: define 

(mXt) := [*{rr,N[0,t]),9>(C;N[0,t])). 

Let IE? q denote expectation in the coupling starting with (r],()- We need also 
to couple the initial configurations. Let JIa be the law of a pair of configurations 
(77, C) with marginal distributions fiA and /1 and such that f](x) = C( x ) f° r all 
x G Z 2 \ Da- It is possible to construct a measure with these properties because 
/i is a product measure. We then have 



l ( , ]A) (i A (m)) -^UU(m)) = / ^(^,C))E ( ,,c)[iAfe)-U((t)] (3.27) 



The number of initial discrepancies is finite, that is, ^ C( x )) < 

|Z?a| < 00. At each site x of Z we have one of three possibilities: (r?(cc) — 
C(x)) = 0, no discrepancies; (?y(a;) — C( x )) + > 0, positive discrepancies; or 
(r](x) — C( x ))~ > 0, negative discrepancies. Following the evolution of the 
particles and the discrepancies we notice that if a positive discrepancy jumps 
over a negative one, then both discrepancies collide, giving place to a coupled 
particle and a hole; if a coupled particle attempts to jump to a discrepancy, the 
jumps occur and then the discrepancy must jump to the site previously occupied 
by the coupled particle. These two behaviors only occur when vertical jumps 
are involved. In the horizontal jumps, discrepancies and coupled particles just 
interchange positions according to the Poisson horizontal (stirring) clocks. 

We say that there is a first class particles at site i at time t when £t(i) = 
r]t(i)Ct(i) = 1, a positive second class particles when (rj^) t (i) = rj t (i) — Ct{i) = 1 
and a negative second class particles when (C 7 ?)* (*) = Ct{i) ~Vt{i) = 1- The first 
class particles occupy initially those sites i occupied by both 77 and (. Locally 
in time, the motion of the first class particles is the one given by generator C v 
but superposed to it there is a pure birth process of first class particles: with 
rate 

P(v()t(h,i2 - 1) (C»?)t(*i)*2) 

the second class particles at (i%,i2 — 1) and (i\ 1 i 2 ) annihilate each other and a 
first class particle appears at (21,12) and an empty site appears at {i\, i 2 — 1). 
Similarly, at rate 

q(v()t(h,i2 + 1) (Cv)t(h,h) 

the second class particles at {i\,i 2 + 1) and («i, i 2 ) annihilate each other and a 
first class particle appears at {i\,i 2 ) and an empty site appears at (ii,i 2 — 1). 

The marginal distribution of a second class particle between two vertical 
jumps (or between a jump and an annihilation) corresponds to the law of a 
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nearest neighbor symmetric random walk — with reflection at the rod when at 
level y — in the horizontal direction. In the vertical direction the motion is 
not Markovian — it depends on the configuration of the first and second class 
particles at the instants of attempted jumps — and either there is an annihilation 
as described above or the second class particles just change horizontal line. For 
instance, at time t, jumps of a (r]() second class particle from site (ii^'a) to site 
(ii,42 + 1) occur with rate 

P(v()t(ii,i2) (1 - €t(h,h + 1)) +q(r)()t(ii,i2)€t(ii,i2 + 1) 

and similarly for the other cases. The first term corresponds to the jump over 
an empty site and the second one to the interchange of positions with a first 
class particle. 

This coupling has the property 

^{E 7 ^) ^ > E 7 M*) ^ = 1. (3-28) 

X X 

i.e., the number of discrepancies cannot increase. 

Since by construction the discrepancies between rj and £ are all located at 
Da, we have the estimate 

%JiA(vt) - MCO] < E E p (^w = z )> ( 3 - 29 ) 

ieD A zeD A 

where X t (t) is the position of a second class particle initially at i. If at site 
i there were no discrepancy we use the convention X\ $ "E 2 (and hence ^ z, 
for all z E Da)- If particles i and j were discrepancies of different sign and 
collided before time t, we also set Xl,X 3 t <£ I 2 . The process X\ has rate 71 to 
move symmetrically in the horizontal direction. If the rod were not present, we 
could dominate JP(X' l {t) = y) by F((X i (i))i = yi), where is the first 

coordinate of the walk. Since without the rod the first coordinate makes just 
a symmetric random walk at rate 71, that probability would be dominated by 

— 1/2 

7 1 times a constant. But with the rod we have to work a bit more. The 
process X z (t) has rate at most 72(p + q) to move in the vertical direction. This 
implies that the time elapsed between the last vertical jump and t is dominated 
by the minimum between an exponential time of rate 72 (p + q) and t. With 
this in hand it is not difficult to prove that also in this case JP(X l (t) = y) is 

— 1/2 

bounded above by 7 1 times a(nother) constant. Here we use that 72 remains 
bounded when 71 goes to infinity. We conclude that for any pair i, y in Da- 

lim JP{X l (t) = y) = (3.30) 

71 too 



Therefore we conclude, combining fl3.25Q , ( fT27| ), J3.29| ) and ( |3.30D and the 
fact that Da is a finite set: 



lim [ dnf y S y (2t)fy = 0. 
71T0O J 



(3.31) 



12 



Remark: We postpone until Section 0, Remark 3, an alternative more 



general proof of Lemma 3.2 which works equally well for a broader class of 
exclusion dynamics (e.g. with speed change) provided the projection of the 
invariant measure on horizontal layers is ergodic for the horizontal dynamics. 

We now prove an intermediate result which is important for the proof of 
Theorem 2.1. 

Proposition 3.1 Let f y be a function depending only on the configuration val- 
ues at the N sites of An(u — 1) or An{]J + 1). Then we have for all t > 0: 

limsupE^ (f Yt (vt) ~ f frM»?idr,)) = 0. (3.32) 

71 Too p \ J / 

Proof: We first want to condition on a sequence T e := (T-f 1 , . . . , T* n ) of marked 
trial jumps before t. Here e G { — 1,+1} is the mark of the jump: +1 for up, 
— 1 for down. Next we consider ai, . . . , a n £ {0, 1} with interpretation on = 1 
if i-th marked trial jump succeeds, ai = if not. Given (T^ 1 , . . . , T^") and 
a := (ax, . . . , «„), we define 

k 

Y^=J2^j. (3.33) 

3=1 

This corresponds to the position of the polymer at time Tt h , given succeeded 
and failed jumps (a%, . . . , ak)- Finally we denote by V p > e := Vjte (rj) the event 
that the polymer in Y p a _ 1 can (for a p = 1) or cannot (for a p = 0) perform the 
jump to Y p c "_ 1 + e p . With this notation, we can write 



n+1 



E (3-34) 

ae{0,l}^- ■••■»> 

X / d^t^i-y^-Tn) (fy- - J /^W^W) • 

Here lJ^\f!f > denotes the monomer distribution at time s = (Tf l n ) + , given the 
successes (ai, . . . , a n ), and IP^o 1 ^ ( a ) denotes the probability of the sequence of 
succeeded and failed jumps prescribed by a at the times T e - The crucial thing 
to realize at this point is that the probability measure ^yZ 2 ^ is absolutely 
continuous with respect to the conditioned Bernoulli measure v p n . In Lemma 



(3.2) below we shall give a uniform bound on the density 



,(71,72) 



du} 

*&' 72) := — tyI - ■ (3.35) 



dv p 



By dominated convergence, the proof of the proposition is reduced to showing 
that for any a 6 {0, l}! 1 ' - -™) and any 5 > 0: 



lim 

71 Too 



d ^yJ )s lL^ 5 ) {frs - J f^vJ" (*»)) = 0- (3.36) 
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The expression inside the limit in the left hand side of (3.36) is bounded by 



lfe 72) lloo pg, 72) W (/ - / f(v)^(d v ))\\ L (3.37) 



Therefore, (3.36) is a consequence of Lemma |3.2| and the following estimate on 



the density *^' 72) . 

Lemma 3.3 Put c(p,x) := [p(x + 1) A p(x — 1) A (1 — p(x))]~ N . For any 
a G {0, 1}^ and for any n £ IN, we have the estimate: 

n-l 

limsup||*£™*>|| < l[c(p,Y p a ) (3.38) 

Proof: We fix a and proceed by induction in n. First put n = 1. By stationarity 
of under the evolution S? \ , we have 

P (.71,72,)' 

^t=v>\Vn- (3-39) 

First consider a\ = 1, i.e., the jump succeeds. Denote V(x) the event that the 
set An(x) contains no monomers. Then we can write: 

f(vy P [d V \vn = J &>T mvm f^v^y ( 3 - 4 °) 



Hence, we conclude 



And we can estimate 



0,(71,72) = i(v(o)) Up ( V ( Yl 1 J1 . (3.41) 



Next consider ol\ = 0, i.e., the jump fails (and thus Y" 1 — 0). We write 

1 



f(vV p [dv\Vn = I v Y ? (d V )[fI(Vr(v))} _ • (3.43) 



Hence, 



So also in that case we have the estimate 

ll*i?' 72 1<^7^<c(p,0). (3.45) 

This proves the claim for n = 1. Suppose the claim is true for n = 1, . . . , p — 1. 
Put a p = 1, the case a p = can be treated analogously. In order to simplify 
the notation, we make some further abbreviations: 
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1 (71,72) ._ 

1. H'a Y" -~ Vp- 



2. v p > 

o ^(71 ,7a) 



(71,72) 



4. 



5 (7i,72) (*) : = sjr™\t) 



r rpCp rpCp-1 

o. lp 1 p -i ■— T P 



6. y° ll >—> a » ~ V „ 



We compute * p 7l ' 72) : 



M/) = (/Xp-r^Vp)) TO 

fdfl p ^S^(T p )(l Vp ) 

fd^ 1 (S^r\r p )(*^) fl Vp 



id** (i v _j 4 7 _ i r ) (^)K 7 -i 72) )) 



(3.46) 



where in the third step we used reversibility of v v p 1 . From ( |3.46 ) we read off 
the density: 



,,(71,72) = H-i S<?F\t p )&™>) 

We first estimate the nominator of the rhs of ( 3.47 ): 
dv p p H-i S^\t p ){^ ) ) 



(3.47) 



A/P- 1 ( c(7l,72)/ W-| \ ,T,(7l,72)\ Vp(Vp-l) 

dvp (Vi (^)(lvp) *p-i J^J" 



> ^p(V^,-i) 

^p(Vp-l) || lT ,(7l,72) 



> 



„ , r , ll^riloollS^^tlv,"^- 1 ^)]!^^-!) 
1 

0(7l), 



<P,Yp-i) 



(3.48) 



where 0(71) tends to zero as 71 f 00 by Lemma 3.2. By the induction hypothesis, 
we obtain from ( gig ), (|0|): 



limsup||^ 7l ' 72) | 

7lT°° 



< limsupH^i^lUcCp,^) 

71 Too 
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p-1 

k=0 



(3.49) 



This finishes the proof of Lemma B.3 and Proposition 3.1 



As a first application we obtain convergence of the one-point marginales of the 



processes {Y t 



(71,72) 



t > 0}. For / : TL — > IR a bounded function, we have, using 



the notation of (2.5). 

(fOQ f(Y ) J* ds (C V J)(Y S 



(3.50) 



By Proposition 3.1 we obtain in the limit 71 f 00: 



Hm e^ f/(y t ) - /(y ) - / d« [(A)^(rfr ? )]/(r s 

lif (/(y t ) - f(Y ) f 

71 Too P \ J Q 



ds(£™f)(Y s ) = 0. (3.51) 



This implies in particular that 



limE^/(^)=Er/« RW ). 



(3.52) 



In order to prove that the processes {Y t 



(71,72) 



t > 0} converge weakly in the 



Skorohod space of trajectories to the random walk {Y^ w : t > 0}, i.e. the 



content of Theorem 2.1, it suffices to show that the process {Y t 



(71,72) 



t > 0} is 



asymptotically Markovian. Indeed, then it is uniquely determined by its single 
time distributions which are those of the random walk {F t RW : t > 0}. More 
precisely it is sufficient to prove the following lemma: 

Lemma 3.4 Let \Tt : t > 0} denote the a-field generated by {{i] s , Y s ) : < s < 
t}. We have 



lim E 

7l Too 



(71,72) 



= 



(3.53) 



Proof: By the Markov property of the process {(rjt, Yt) : t > 0}, 



E 



(f(Xt-,)) 



f(Y s )+TE 



(71,72) 

x 6y„ 



C v J(Y r )dr. (3.54) 



Therefore, it suffices to show that 

r t-s 



lim E 

71 Too 



(71,72) 

7] s XSy s 







£„ r /(y r )eZr- / dr / £ v f(Y r )vf(dr))\ =0 (3.55) 



Since the trial jumps of the polymer are on the event times of a Poisson process 
with rate independent of (71,72), we can write 



C n J(Y r )dr = 



r+e 



£ Vrl f(Y r ,)dr' +o(e), 



(3.56) 
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where o(e) goes to zero in L 2 (JP^ 1 ^ 2 ), uniformly in (71,72), when e tends to 
zero. Therefore, it is sufHcient to show that 

fy r (Vr)dr- / rf s {drj)f Ys (v) 



H™ I U? (7i .72) lei (71. 72) 



(3.57) 



Following the same strategy as in the proof of Proposition [0|, i.e., by estimates 
on the density of the monomer distribution with respect to the appropriate 
conditioned Bernoulli measure, this reduces to showing that for any e > 0, for 
any y G Z and for f y depending on layer y + 1 or y — 1: 



lim IE 

7! too 



(71,72),!/ 



dsfy(Vs)- / v y P (dv)fy(v) =0 



(3.58) 



Putting /y := /j, — the expression inside the limit in the left hand side of 

( [3.58 ) can be rewritten as 



V v(drj)- 



ds 



drf y S v hirf2) (\r-s\)f y 



< 



~ ds dr\\f y \\ L 4Sl^ 2) {\r-s\)f v \\ 



L 2 ■ 



(3.59) 



Hence we obtain ( p. 58 ) as an application of Lemma |37 
Arrived at this point, we know that any weak limit point of the processes 
1^(71,72) ■ t > 0} equals in distribution the random walk {Y t RW : t > 0}. 
Hence, to finish the proof of Theorem 2.1, it is sufficient to see that such a weak 
limit point actually exists. This is an easy task: 



Lemma 3.5 The sequence of processes {Y t 



(71,72) 



t £ [0, T], }( 71l72 ) is tight. 



Proof: Since the number of jumps the polymer makes in [0,T] is bounded by a 
mean one Poisson process, we have 



and also 



2T 

P( sup |Y s (7i.72)| > M) < — , 
0<t<T M 



HmP( SUp |y s (7i>72) _ F( (7i,72)| > e) = 

fi J-° s,te[0,T],\s-t\<6 



(3.60) 
(3.61) 



This proves tightness (cf. Theorem 1.3 p. 51 of 0). 



4 Additional remarks 



Remark 1: What happens when the system is out of equilibrium? For instance, 
start the monomers in a homogeneous product measure. When the density is 
constant and equal to p G [0, 1] (no p,q,i2— dependence in (2.6), the measure 
v v is no longer invariant for the monomer dynamics S v (t) at fixed rod position 
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y. However in the limit 71 — » 00 the polymer will perform a continuous time 
random walk with rates a(l — p) N and 6(1 — p) N for up and down jumps re- 
spectively. Significant corrections in the case 71 < 00 can be expected, cf. 0. 
Another problem is obtained if we start the monomers from a sharp density pro- 
file. That is, above the polymer the fluid density is constant p\ and under the 
polymer the density is also homogeneous equal to p2- In this case the vertical 
density will follow a discrete space noiseless Burgers equation: p(i,t) e [0,1], 
t e H, i £ 1 

= - pp (i,t)(l-p(i + l,t)) - qp(i,t)(l-p(i-l,t)) 

+ pp(i-l,t)(l-p(i,t)) + qp(i+l,t)(l-p(i,t)) 

with initial condition p(i,0) = pil(i < 0) + p\I{i > 0). The limiting motion of 
the rod will be a non-homogeneous (in time) Markov process described by 

= a[l-p(Y t + l,t)] N [f(Y t + l)-f(Y t )} 

+ b[l-p(Y t -l,t)] N [f(Y t -l)-f(Y t )]. 

where T t is the sigma field generated by {Y s : s < t}. These results can be 



obtained with the techniques we used to prove Theorem 2.1 and will be the 
content of a future publication, cf. || . 

Remark 2: One may wonder how general the results are. As an illustration of 



this we consider the following somewhat abstract modification of Lemma 3.2 
Suppose that p is a reversible measure on {0, 1} both for a monomer dynamics 
with generator C\ and one with generator £2 . As an example, we could keep in 
mind the case where £12 = £1 + £2 is a Kawasaki dynamics (exclusion process 
with speed change) at finite temperature with C\ generating the horizontal 
and £2 generating the vertical jumps; p is the corresponding Gibbs measure. 
The measure p is then also reversible for Cj 2 = 7£i + £2- Now we insert the 
polymer and we fix it at some position y £ TL. The dynamics of the monomers 
is now conditioned on having no monomers in the excluded volume A^{y): 
t]t(i) = 0, Vi £ Ajv(y),Vt > 0. The new generator is £j / 2 7 = + £2 obtained 
by setting all of the original rates equal to zero for all updating that would cr eate 
a mo nom er in the region A^{y) (the dir ect analogue of what was done in ( |2.l|) 
and (f2~2f)). It follows then from Lemma |1] that p v = p(-\rj(i) — 0, V« £ A N (y)) 
is reversible for C\ 2 . We finally denote by p v x the restriction of p v to the layer 
at height x (i.e., the set {i S Z 2 ,«2 = %})■ This measure is reversible for C\. 
We have the following result: 

Proposition 4.1 Denote by S^it) the semigroup with generator C\£ ■ Assume 
that for all x p v x is ergodic for L\. Let f x be a function in L 2 (p y ) with depen- 
dence set on layer x ^ y. We have: 

km \\Sq(t)f x - [ dpi = 0. (4.62) 

7T00 " J 

Proof: By ergodicity C\ has simple eigenvalue with corresponding eigenspace 
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the constant functions. Hence by the spectral theorem, 



(4.63) 



where IE^ 1 . denotes the spectral measure of the selfadjoint operator C \ . There- 
fore, we have to show that if f x is a function on layer x such that 



then 



lim||SJ(*)|| i2(MH) =0. 

7*00 



(4.64) 
(4.65) 



For every tp in the domain of 1 



lim -(C\£ip) = C[ip. 

7T00 7 



(4.66) 



Hence the spectral measures JEj ~"j 12 converges weakly to the spectral measure 



IE^. Therefore, we can estimate 



\\s^{t)f x \\l w 



< 



< 



— A. /*f '7 
Jx ,Jx 



(dX) 



,—yts 



(4.67) 



Letting 7 tend to infinity, and then S to zero, using (4.64), we obtain (4.65). 



Remark 3: Proposition 4.1 is general but has a strong hypothesis: the ergod- 
icity of the one-layer horizontal dynamics. This is known only in a few cases, 
in p articular in the symmetric simple exclusion process we treated in Lemma 
3.2. It is however expected to be true at least for high temperature Kawasaki 
dynamics. 
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